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In a previous paper we proposed a Projected Configuration Interaction method that uses sets of 
axially deformed single particle states to build up the many body basis. We show that the choice 
of the basis set is essential for the efficiency of the method, and we propose a newly improved 
algorithm of selecting the projected basis states. We also extend our method to model spaces that 
can accomodate both parities, and can include odd-multipole terms in the effective interaction, such 
as the octupole contributions. Examples of *^Fe, ^''Ni, ®*Se, ™Se and ^®Se are calcualted showing 
good agreement with the full Configuration Interaction results. 
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I. INTRODUCTION 

The full configuration interaction (CI) method [U, 
using a spherical single particle (s.p.) basis and realis- 
tic Hamiltonians, also know as the nuclear shell model, 
has been very successful in describing various properties 
of the low-lying states in light and medium nuclei. The 
main limitations of this method are the exploding dimen- 
sions with the increase of the number of valence nucleons, 
or/and with the increase of the valence space. Although, 
there are continuous improvements to the CI codes [3, 13| 
and computational resources, the exploding CI dimen- 
sions significantly restrict the ability to investigate heavy 
nuclei, especially those which exhibit strong collectivity. 
The deformed mean-field approaches, however, have the 
ability to incorporate the collective effects at the single 
particle level. The mean-field description in the intrinsic 
frame naturally takes advantage of the spontaneous sym- 
metry breaking. This approach provides some physical 
insight, but the loss of good angular momentum of the 
mean-field wave functions makes the comparison with the 
experimental data difficult. The CI calculations in spher- 
ical basis provide the description in the laboratory frame. 
The angular momentum is conserved, but the physical in- 
sight associated with the existence of an intrinsic state 
is lost. One important aspect of the CI approach is its 
ability of using all components of effective interactions 
compatible with a given symmetry, but restricted to a 
chosen valence space. Examples of realistic Hamiltoni- 
ans, such as the USD [^in the sd shell, the KB3 [6j, 
FPD6 and GXPFl Q in the pf shell, have provided a 
very good base to study various nuclear structure prob- 
lems microscopically. 

The recent history of projection techniques combined 
with CI particle-hole configurations includes the Pro- 
jected Shell Model(PSM) [9, lOj and the Deformed Shell 
Model (DSM) proposed in reference PSM uses a de- 
formed intrinsic Nilsson-|-BCS basis projected onto good 
angular momentum, and a multipole-multipolc Hamilto- 
nian that diagonalized in the space spanned by the pro- 
jected states. The Nilsson model [I4I has been proven 



to be very successful in describing the deformed intrin- 
sic single particle states, and the quadrupole force was 
found to be essential for describing the rotational mo- 
tion 13]. PSM was proven to be a very efficient method 
in analyzing the phenomena associated with the rota- 
tional states, especially the high spin states, not only for 
axial quadrupole deformation, but also for the octupole 
and triaxial shapes [ij, However, its predic- 

tive power is limited because the mulitpole-multipole plus 
pairing Hamiltonian has to be tuned to a specific class of 
states, rather than an region of the nuclear chart. The re- 
cently proposed DSM is using the same realistic effective 
Hamiltonian as the full CI method, and a Hartree-Fock 
procedure to select the deformed basis. One can only as- 
sume that this procedure would not be very accurate for 
quasi-spherical nuclei. The main difficulties for all these 
models is a proper selection of the deformed basis. Their 
accuracy can only be assessed by comparison with the ex- 
act results provided by the full CI method using the same 
effective Hamiltonian, and not by direct comparison with 
the experimental data. Other model using similar tech- 
niques includes MONSTER, the family of VAMPIRs [ll, 
and the Quantum Monte Carlo Diagonalization (QMCD) 
method [l7j . 

In a previous paper [Tsj we proposed a new method 
of calculating the low-lying states in heavy nuclei using 
many particle-hole configurations of spin-projected Slater 
determinants built on multiple sets of deformed single 
particle orbitals. This Projected Configuration Interac- 
tion (PCI) method takes advantage of inherent mixing 
induced by the projected Slater determinants of varying 
deformations with the many particle-hole mixing typical 
for the Configuration Interaction (CI) techniques. Di- 
rect comaparison between PCI and CI results are always 
possible, provided that the deformed s.p. states are al- 
ways obtained starting from a valence space of shperi- 
cal orbitals. In Ref. we use a simple mechanism 

of selecting a number of basic deformed Slater determi- 
nants in the sd and pf model space, denoted \Kj, >, 
by searching for the minimum energy of fixed configu- 
ratin of deformed s.p. orbitals. Starting from each basic 
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deformed Slater determinant, a number of particle-hole 
excited configurations were considered under some selec- 
tion criteria (see Eq. (23) of Ref. [ll]) to keep the total 
numebr of basis states manageable. Having the deformed 
basis of Slater determinants chosen, one can use standard 
spin projection techniques to solve the associated eigen- 
values problem. The method proved to be very accurate 
in Ohio model spaces, such as sd and pf where one can 
easily keep track of different deformed orbitals. The dif- 
ficulties usually appear for quasi-spherical nuclei, such as 
^^Ni, when special attentions has to be given to the se- 
lection of the basics states \Kj, >. A simlar problem 
arises for the case of mixed parity valence space, such 
f5pg9 (see below), due to difficulties in tracking fixed 
configurations of nucleons filling the s.p. orbitals around 
the level-crossing deformations. 

The paper is organized as follow. Section II presents a 
brief outline of the PCI formalism that was expanded in 
Ref. [18]. The new algorithm to select the PCI basis is 
discussed in section III. Section IV analyzes the efficiency 
of the new method in the case of the quasi-spherical nu- 
cleus ^^Ni. Section V is devoted to the study of several 
nuclei that can be described using the mixed parity va- 
lence space f5pg9. Conclusions and outlook are given in 
section VI. 



determined by the quadrupole £2, hexadecupole £4 as in 
Ref. [Ill, but also octupole £4, etc. 

The general form of the nuclear wave function is taken 
as a linear combination of the projected SDs (PSDs), 

mAl)=J2flK.HlK\^)> (4) 
Kk 

where Pljj^ is the angular momentum projection opera- 
tor. The energies and the wave functions [given in terms 
of the coefficients ffj^,^ in Eq.lU)] are obtained by solving 
the following eigenvalue equation: 

{Hkk.K'h' ~ EfNj^^ji,^,)ffj(^, = 0, (5) 

K'k' 

where Hj^^ and j^,^, are the matrix elements of 
the Hamiltonian and of the norm, respectively 

hl,k'.' = mp'KK'W), (6) 

N'kk.k:- - {'^\PkK'W)- (7) 
More details about the formalism can be found in Ref. 




II. THE METHOD OF THE PROJECTED 
CONFIGURATION INTERACTION (PCI) 

The model Hamiltonian used in CI calculations can be 
written as: 



H 



E 



ClCk, 



(1) 
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where, c\ and ci are creation and annihilation operators 
of the spherical harmonic oscillator, and Vijki are one- 
body and two-body matrix elements that can be obtained 
from effective interaction theory, such as G-Matrix plus 
core polarization ^19*1 , which can be further refined using 
the experimental data 8, 20]. 

One can introduce the deformed single particle (s.p.) 
basis, which can be obtained from a constraint Hartree- 
Fock (HF) solution, or from the Nilsson s.p. Hamiltonian 
The deformed s.p. creation operator is given by the 
following transformation: 



(2) 



where the matrix elements Wki — {bk\ci) are real in our 
calculation. The Slater Determinant (SD) built with the 
deformed single particle states is given by 



(3) 



where s refers to the Nilsson configuration, indicating the 
pattern of the occupied orbits, and e is the deformation 



III. CHOICE OF THE PCI BASIS 

The analysis made in Ref. ,1^ indicated that one of the 
most important problem of the PCI method is a proper 
selection of the PCI basis. As introduced in our previous 
work [1^, the general structure of the PCI basis is 



'Op — Oh, np — nh 

|ki,0), \ni,3),--- 

|k2,0), 1^2, j),- •• 

|KAr,0), |KAr,j),--- 



(8) 



where [k^, 0) (i = 1, ...N) is a set of starting states of dif- 
ferent deformations. Assuming that we've found these 
|k, 0) SDs (skipping the subscript i to keep notations 
short), relative np-nh SDs, |K,j), on top of each |k, 0) 
are selected using the constraint [18] 



(Eo - EjY 



4|y|2) >K 



(9) 



where E^ = (k, OjiJjK, 0), E.j = (k, jJiJjK, j), V = 
(k, 0]i?|K,j) and Ecut is a parameter. 

The Jk, 0) SDs need to be properly chosen in order 
to get good accuracy. In our previous work (l8| . we 
have chosen the SDs with the lowest unprojected expec- 
tation energy for each configuration, and we used the 
same basis for all the spins. That approach proved to 
work well for quite deformed nuclei, limiting its range 
of application. For instance, the description of ^^Ni with 
GXPFIA [iO] exhibits a spherical ground state minimum 
which is selected as a basis SD. This spherical SD has a 



3 



good spin 7 = 0, which won't be useful if we calculate 
/ ^ states. This example suggests that a more effi- 
cient method may involve choosing different basis sets 
for different spins. Another problem with this method of 
selecting the |k, 0) basis states is that for each configu- 
ration, only one SD was selected. Therefore some states, 
such as the /3— vibrational states, cannot be described 
unless two or more shapes for the same configuration are 
artificially included outside of any algorithm. 

To address these problems, we developed a new 
method of finding an cfhcient set of |k, 0) states. The 
deformed single particle states are generated from the 
Nilsson Hamiltonian shown in Eq. (2) of Ref. [T^. For 
simplicity, we set 



(10) 



In a first step, at each deformation e = (62,63, €4), we 
build many Slater determinants(SDs) denoted by |s,e), 
where s denotes a configuration of nucleons occupying 
the deformed single particle orbitals. These SDs are pro- 
jected onto good angular momentum /, and the projected 
energy is calculated 



KK\ 



(11) 



We then identify the configuration Sa which has the 
lowest Eaxpil , s, e) at each shape e. Searching over all 
possible deformations e, we obtain the energy surface of 
EcxpiI,Sa,e) as a function of e. The SD \sa,ea) which 
has the lowest Ecxp{I,Sa,e) is chosen as the first \k,0) 
state denoted as 



|k1,0) = \Sa,ea). 



(12) 



The next step is to find the second |k, 0) state. We try 
all possible |s,e), and for each |s,e), we build the 2x2 
Matrix pair (A, B), 



^ ^ I'i/n Hi2\ g ^ (Nn iVi2 

H21 H22 I ' \N21 N22 



where 

H,, = {i\HPi,i,\j),N,, = {ilPiiKlj), (13) 
with |z(j) = 1) = 1^1,0), |*(j) = 2) = |s,e). (14) 

Solving the generalized eigenvalue problem 

Aa; = XBx, (15) 

we get two eigenvalues, Ai and A2, and their sum, 

^2 = Af^+Af. (16) 

The SD jsf), £{,) with the lowest 52 is selected as the second 
\k, 0) denoted as 



|k2,0) = |sfc,efc). 



The process of finding more |k, 0) SDs can be continued 
in a similar manner. Assuming that we have found the 
(n — 1)— th |k, 0) SD, |k„_i,0), then |k„,0) is chosen as 
the \sx,£x), corresponding to the lowest 5'„. Here, 
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Ai"). 



(18) 



and a("\a("\..., Ar are eigenvalues of Eq. with 
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HPirK\j),N,, = { 




(19) 



K(j)) = |Ki(j),0), if 



1,2, 



(20) 



Sometimes we may only use part of the Sn sum over 
the lowest A^'s as a selection criteria, i.e.. 



5fc = A(")+A(")-f ••• + Ai"\(l<fc<n), 



(21) 



and = Sn by definition. 

Evaluating the projected energies for all SDs takes a 
long time to calculate, and not all of them may be neces- 
sary. Therefore, we enforce additional truncations. First, 
the HE energy, i?HF, is calculated. Next, at each shape e, 
the SDs having all particles occupying the lowest Nilsson 
orbits are considered as the Op-Oh SDs. All particle-hole 
excitations up to 4p-4h built on these Op-Oh SDs are cre- 
ated, and their expectational energies 



E, 



'cxp('5i ^) ^ 



e\H\s, 



(22) 



are calculated. Those SDs satisfying i?cxp(s,e) — E^p < 
where i?cxpup is a input parameter, are saved. Ei- 



cxpup 1 



(17) 



E, 

nally, the projected energies Ecxp{I, s, e) of the saved SDs 
are evaluated, and compared with the lowest projected 
energy Eexp{I,Sa,ea) available. We keep those SDs sat- 
isfying Eexp{I,s,e) - Eexp{I,Sa,ea) < £^pjup(/), whcrc 
i?pjup(/) is a input parameter, and we discard the oth- 
ers. The values of the parameters Scxpup and ii^pjup(/) 
must be large enough so that the |k, 0) can be properly 
chosen, but too large i?cxpup and i?pjup(/) value may re- 
sults in wasted computation without any improvement in 
accuracy. 

Here we summarize the advantages of the new method 
of selecting the basis states \k,0). Firstly, as already 
mentioned, the method proposed in Ref. [1^ uses the 
same |k,0) states for all spins, however, certain |k,0) 
states may not bring any contribution to certain spins. 
The new method improves the efficiency of the PCI basis, 
by choosing different \k,0) states for different spins. 

Secondly, the \k, 0) chosen by the new method may in- 
clude two or more shapes for the same configuration s. 
Therefore, the present method explicitly includes the idea 
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imbedded in the Generator Coordinate Method (GCM) 
PH . and may be used to describe some collective vibra- 
tions, such as the /3— vibration. 

Thirdly, there are no limitations on the K values. In 
Ref. we only selected basis states with relatively 

small- if values to keep the basis dimensions manageable. 
This limitation could be a problem in the case of the high 
spin states, for which the high-if configurations could 
be close to the yrast line. The new method described 
here selects basis states with all possible SDs satisfying 
\K\ < I. 

Finally, and perhaps more importantly, the drawback 
of getting large overlaps between different basis SDs, 
which is typical for an uncorrelated selections of the basis, 
is avoided by construction in the new method. Getting 
basis states with large overlaps leads to many spurious 
states due to the zero eigenmodes of the norm matrix. 
The significance of these zero modes is that some of the 
basis states that exhibit large overlaps bring insignifi- 
cant contribution to the solutions of the Hill- Wheeler 
Eq. (5), while unnecessarily increasing the dimensions 
of the problem. Those useless SDs can be automatically 
filtered out by the present method because the overlap 
problem has been fully considered step-by-step when the 
generalized eigenvalue equations are solved. 



IV. CALCULATIONS OF *NI 

Using the new method, we recalculated the nuclei ^®Ni 
with GXPFlA interaction^]. Let's first consider the 
case of / = 0. Both €2 and €4 span the interval from —0.45 
to 0.45 in steps of 0.03. The first basis state, 0), hav- 
ing the lowest projected energy, corresponds to / = 0. In 
Fig. [H the surfaces of the _Boxp(i, Sa, e) (See Eq. (fTTj) ') 
and i?cxp(sa,e) (See Eq. ([22]) ) are plotted as functions 
of €2 and €4. Our calculation shows that the configura- 
tion Sa has all 16 valence particles in ^^Ni occupying the 
orbits coming from the I/7/2 subshell. The unprojected 
minimum. Fig. (l)(a) is at —203.800 MeV, and its shape 
is spherical, consistent with the HF result. However, the 
projected energy surface presents a quite different pic- 
ture. There are four minima around the spherical shape, 
the lowest one has —204.473 MeV and a small oblate de- 
formation, £2 — —0.09 and £4 = —0.09. This energy is 
673 keV lower than the HF energy, and 1.236 MeV above 
the exact CI ground state energy of —205.709 MeV. This 
K,0} state is a good candidate for the ground state. 

The second basis state, |k2, 0), has the same configura- 
tion as the first one, but a different shape characterized 
by £2 — —0.24 and 64 = —0.15. Comparing with Fig. [U 
this shape is quite different from any of the remaining 3 
minima. The reason is that the SDs at those 4 minima 
are highly overlapping each other after the angular mo- 
mentum projection. Once the lowest one is picked up, the 
others will automatically be filtered out by the present 
method. The |k2,0) corresponds to the first excited 0^ 
state, which might be called a /3— vibrational state. 
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FIG. 1: (Color online) Unprojected energy surface (left panel) 
and the projected energy surface with 7 = (right panel) for 
the ground state of ^^Ni with the GXPFlA interaction. The 
lowest energy is marked by '©' 
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FIG. 2: (Color online) Si - Si-i values of 0) at / = for 
^®Ni as a function of £2. £4 are included in the calculation. 
|aii,0) is the lowest one. 



The third basis state, |k3, 0), has a prolate shape with 
£2 = 0.27 and £4 = 0.06. Its configuration can be ob- 
tained starting from 0), but with 4 particles jumping 
from the \n\ = 7/2 (I/7/2) orbits to the \n\ = 1/2 (2^3/2) 
orbits. |k3,0) can generate a deformed rotational band, 
which has been observed in experiments ^23|. Here, we 
only create the band head, which is the third 0+ state. 
Informations about higher |k, 0) SDs are shown in Fig. 
[2] The value Si — Si-i indicates the energy position of 
each |Ki,0) state. 

Once we have selected the \k,0) SDs, we perform the 
PCI calculations. There are two parameters used in PCI: 
one is the number of |k, 0) SDs, n, and the other is the 
i?cut used in Eq. (9) to select the number of particle-hole 
excitations on top of each |k,0) [Tsj . 

It is interesting to study how many |k,0) are needed 
to describe the low- lying states. Fig. (3) (a) shows the 
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FIG. 3: (Color online) Left panel, 7 = PCI energies for '^^Ni 
as functions of n with E^ut = 1 keV. Right panel, 7 = PCI 
energies for ^^Ni as functions of 7Jcut = 1 with n = 8. Full CI 
results are also shown as the open circles by comparison. 




FIG. 4: (Color online)The lowest 5 energies at each spin for 
^®Ni calculated using PCI (open circles) and full CI (filled 
circles) . 



PCI energies as functions of n for Ecut — 1 keV. For 
n — 1 the PCI dimension is only 180, yet the first O"*" 
PCI energy is -205.409 MeV, only 300 keV above the 
exact full CI value. |k2,0) must be included to describe 
the second 0"*" state, and |k3,0) reproduces the third 0+ 
state. To accurately describe more excited states more 
\k, 0) SDs are need. For the lowest 5 states in ^^Ni a good 
approximation can been achieved starting with n = 7. 
If one keeps on increasing n, then those 5 energies will 
become closer and closer to the CI values. For instance, 
with n = 15, the the ground state PCI energy becomes 
-205.603 MeV, just 100 keV above the exact value. 

The PCI energies are also affected by the i?cut param- 
eter. For Ecut — 1000 keV no particle-hole excited SDs 
are included, and the PCI dimension is the same as the 
number of |k, 0) SDs included, n = 8 in this case. There- 
fore the PCI energies are exactly the values of Xi in Eq. 
P^ . More particle-hole excited SDs can be included by 
reducing the value of i?cut- For example, by decreasing 
from £;cut = 1000 keV to E^nt = 1 keV, the PCI energies 
drop ~ 1.0-^2.3 MeV for the lowest 5 states, and become 
close to the full CI values. By decreasing from Ecut — 1 
keV to Ecut — 0.2 keV the energy drop becomes slower, 
and is around ^ 100 200 keV. The PCI energies for 
/ 7^ states are also calculated, and shown in Fig. |4l the 
number of |k, 0) is n = 15 and Ecut = 1 keV. One can 
observe good agreements between the PCI and the CI re- 
sults, including the states in the rotational band starting 
at about 5 MeV. 



V. CALCULATIONS IN f5pg9 VALENCE SPACE 

We have also extended our calculations to the f5pg9 
valence space, which includes the I/5/2, 2^3/2, 2pi/2, and 
the lgg/2 spherical shells. 159/2 orbital has positive par- 
ity, and the other fp orbitals have negative parity. There- 
fore, SDs with both parities can be built for any num- 



ber of nucleons. The positive parity SDs are those with 
even number particles occupying the fp orbitals, while 
the negative parity SDs have odd number particles oc- 
cupy the fp orbitals. The angular momentum projection 
does not change the parity. The parity of the projected 
states remain the same as that of the original SDs. There- 
fore, one can split the PCI basis into the positive parity 
part, and the negative parity part, /~, at each spin 
/. With the present method, different |k, 0) SDs can be 
generated separately for the /+ basis and I~ basis. 

The interaction for f5pg9 shell space was taken from 
Ref. [2^. It includes, besides the usual qqudrupole, hex- 
adecupole, and pairing terms, octupole contributions and 
monopole corrections. In all cases, 20 \k, 0} SDs are taken 
for each basis and E'cut was fixed to 1 keV. The first 
example we analyze is the N=Z nucleus ^*Se, which is 
known to be deformed with competing oblate and pro- 
late deformations. The energies of the 20 \k,0) SDs for 
both = 0+ basis and = 0^ basis are shown in Fig. 

m 

It is known that ^*Se nucleus exhibits shape coexis- 
tence features. The constrained HF calculations of Ref. 
[2^ as well as our results in Fig. [Hlshow that there are two 
minima. Both of them are axially and reflection symmet- 
ric. The lowest minimum has —40.718 MeV and oblate 
shape, and the second one has —39.956 MeV and pro- 
late shape. It is interesting that the results of our new 
method presents the same picture as one can observe in 
the left panel of Fig. [5l where the lowest \ki,0) is oblate 
with 51 = Ec^p{0,Sa,ea) = -42.405 MeV and |k2,0} is 
prolate with 52 — 51 = —41.549 MeV. However, both en- 
ergies are about 1.6 MeV lower than those obtained by an 
HF procedure, due to the angular momentum projection. 

For the 0~ basis the lowest |k, 0) SD lies at 
— 36.44MeV, which is relatively high (see FigO, and has 
a prolate shape. The corresponding configuration is the 
same as that of the second HF minimum, except that 



6 




-0.4 -0.2 0.0 0.2 0.4 



-0.4 -0.2 0.0 0.2 0.4 




FIG. 5: (Color online) Si - Si-i values of \Ki,0) for r = 0+ 
(left) and — 0~ (right) in ^*Se as a function of £2- £4 are 
included in the calculation. 



FIG. 7: (Color online) The lowest 3 energies at each 
spin/parity for ®*Se calculated using PCI (open circles) and 
full CI (filled circles). 



CHF for ^^Se(MeV) 




TABLE I: PCI Dimensions compared with those of full CI for 

™Se 
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FIG. 6: (Color online) Energy surface provided by Con- 
strained Hartree-Fock (CHF) calculations as a function of 

92 = ^ (i)^ Y20 and qs = {^)^Y3o. (6 is the harmonic 
oscilater length) 



one particle was excited from the f2 = 8/2(^3/2) orbital 
to the n = 3/2(39/2) orbital. Because ^^Se has N ^ Z, 
this excited particle can be either a neutron or a proton, 
and therefore there are two different \k, 0) SDs with the 
same shape and the same low energy. Similar cases ap- 
pear for other |k, 0) SDs. Therefore, in Fig. \5lb) each 
symbol corresponds to two different |k, 0) SDs. The po- 
sition of the second lowest symbol is only about 130 keV 
above the lowest one, and its configuration is similar to 
the lowest one, but with the odd particle excited from 
the f2 = 1/2(^3/2) orbital to the r2 = 1/2 ((79/2) orbital. 

Using the |k, 0) states described in Fig. [5l one can 
calculate the PCI energies for the 0+ and 0~ states, which 
turned out to be very close to the CI results. For 1^0, 
similar good results have also been achieved, as shown in 
Fig.m 

As indicated in the above discussions, the PCI method 



not only provides a good approximation for the CI re- 
sults, but it is also a convenient tool to gain some insight 
into the physics of the nuclear states. One interesting 
example is the lowest state in Fig. El^b), which is a 3^ 
state. As shown in Fig. [51 the (lowest) |ki,0) SD has 
K'^ — 3^ and oblate deformation. The configuration of 
this 0) is the same as the oblate HF minimum, except 
that one particle was excited from the H. — 3/2(^3/2) or- 
bital to the ri = 9/2((79/2) orbital to form a = 3~ 
SD. The second SD, \k2, 0), has the same energy and the 
same shape as |ki,0) because N = Z, and due to the 
isospin symmetry of the adopted Hamiltonian. If only 
the particle- hole excitations built on |ki,0) are included, 
one obtains a PCI energy of -40.469 MeV. If the |k2, 0) 
SD is further included, the PCI energy drops 300 keV 
to -40.769 MeV, This energy is only 300 keV above the 
exact value of —41.043 MeV. However, the PCI energy 
for 71 = 20 is -40.843 MeV, only 70 keV lower than what 
one can obtain with n = 2. Therefore, it is clear that the 
lowest 3~ state has mostly contributions from the lowest 
2 oblate K'^ = 3" SDs, i.e. |ki,0) and |k2,0). 

We have also calculated states of both parities in ™Se. 
The results are shown in Fig. [HI Once again, the PCI 
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FIG. 8: (Color online) (a), Si-S,-i values of \k, 0> at r = 3" 
in ^*Se. (b), PCI energies (Open Symbols) and CI energies 
(Fill Circles) at r = 3" in ''^^Se. The open Symbols from left 
to right refer to PCI calculations with n = 1, 2 and 20. 
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FIG. 9: (Color online) The same as Fig. [7] but for ™Se. 



results are very close to those of full CI for both positive 
parity and negative parity for a wide range of spin values. 
Hovifever, much smaller dimensions of the PCI matrices 
are necessary. The PCI dimensions correponding to the 
/■^ calculations in Fig. [S] are compared in Table |T] with 
the full coupled-I CI dimensions. The PCI dimensions are 
small fractions, roughly 10"'', of the full CI dimensions. 
As is well known, the most serious problem with full CI 
method is the explosion of the dimensions as the number 
of the single-particle valence states, and/or number of 
valence nucleons. However, this problem seems to be less 
of an issue for the PCI method. The total PCI dimension 
can be estimated as the product of two numbers, n x m, 
where n is the number of |k, 0) states and m is the number 
of particle- hole excitations selected by Scut in Eq. 
Our investigations indicate that n is related to how many 
low-lying states of a given spin one wants to accurately 
describe. For instance, if we are interested in only the 
yrast state, quite often a good approximations can be 
obtained with n = 1 or 2. n — 20 seems to be enough to 
describe the lowest 3-5 states of each /'^ in the present 



FIG. 10: (Color online) Lowest 3 0+ energies of ^'^Ge and '^^Se 
calculated by PCI (filled symbols) and full CI (open circles) 



calculations. As regarding m, the \K,j) SDs are limited 
to Ip-lh and 2p-2h excitations according to Eq. ([9]). 
Note that \K,j) has the same K'^ as that of \k,0). As 
the parameter i?cut is enforced via Eq. ([9]), m can be 
significantly reduced. For instance, in the case of I'^ = 0+ 
for ™Se, m2p2h = 801 for |ki,0), but only 138 SDs were 
finally included if Ecut — 1 keV. 

The exploding CI dimensions has as a consequence a 
rapid increase of the computing time necessary for full CI 
calculation. Using the modern coupled-I code NuShellX 
4] , the full CI calculation of all states in Fig. [5] could 
take almost one year when only one processor is used. 
The calculation of the lowest 3 states of each J'^ in ''°Se 
would take in average about 20 days. For the same calcu- 
lation, PCI takes around 5 hours for each I^. The main 
computational workload in PCI is related to the calcula- 
tion of the dense matrices H, and N in Eq. It should 
be metioned that extra time is needed to to search for the 
optimized set of |k, 0) SDs. The computing time can be 
affected by: (1) the number of mesh points used for the 
shape paramters; (2) the values of the parameters i?oxpup 
and i5pjup(/'^) that decides how many SDs are considered 
in the optimization process; (3) the total number n of 
|k, 0) basis states selected; (4) the value of -Ecut- For ex- 
ample, in the calculation of /'^ ~ 0"*" in ''"Se, both €2 and 
£4 run from —0.3 to 0.3 in steps of 0.02, -Boxpup = 7 MeV 
and £;pjup(/'' = 0+) = 5 MeV, and E^cut = 1- Under 
these conditions, it will takes about 10 hours to obtain 
20 \k, 0) SDs using one processor. For other I"^ , the com- 
putational time ranges from few hours to 1 or 2 days. 
However, the total time for a PCI calculations is at least 
10 times shorter than that of the corresponding full CI 
calculation for the case of ™Se. 

Finally, we used the new PCI method to calculate the 
low lying O"*" states in ''^Ge and ''^Se. The nuclear struc- 
ture of these two nuclei is relevant for the double beta 
decay (DBD) process of '^^Ge. DBD is one of the most 
actively investigated nuclear physics problem, which may 
reveal new physics beyond the Standard Model, includ- 
ing the absolute values of the neutrino masses. Full CI 
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calculations [U [2^, of the 2-neutrino and neutrino- 
less DBD matrix elements have been carried out for some 
DBD nuclei up to ^^^Xe. However, for heavier DBD nu- 
clei ^^°Nd and ^^^U, the huge CI dimensions make the full 
CI calculation unmanageable. PCI can take full advan- 
tage of the deformation, and an efficient truncation could 
be obtained for well deformed nuclei, such as ^^°Nd and 
^^^U. As a first inroad into this problem, the low-lying 
0+ states ^^Ge and ""^Se are calculated using the present 
version of the PCI, and are compared with full CI re- 
sults in Fig. [ini Using only 6 |k, 0) SDs {n — 6) for 
each nucleus, the PCI dimensions are 561 and 647 for 
''^Ge and ''^Se, respectively. The calculated PCI energy 
of the lowest 0+ state for ^^Se is 200 keV higher than the 
exact value, and only 86 keV higher for ^^Ge. In addi- 
tion, good approximations for the excited 0"'' states have 
also been reached. Given these encouraging results, one 
would hope that PCI calculations could be successfully 
performed for the heavy deformed DBD nuclei, such as 
^^°Nd and ^■^'^U, in a not so distant future. 

VI. CONCLUSIONS AND OUTLOOK 

In this article we propose a newly improved algorithm 
of selecting the basis of Slater determinants that can 
be used with the Projected Configuration Interaction 
method introduced in Ref. [ll]. The new algorithm de- 
pends on a number of parameters that can be used to fine 
tune its efhciency. Its main advantages over the original 
method of selecting of the basis are summarized at the 
end of Section III. 

We used the new algorithm to revisit the calculation 
of ^^Ni, quasi-shpcrical nucleus that has a relatively low- 



lying rotational band. We were able to calculate its low- 
lying states very efficiently, and with good accuracy, gain- 
ing also insight into the physics of these states. We have 
also use the new method to analyze some Se and Ge 
isotopes in the f5pg9 valence space. Both natural and 
unnatural parities can be accurately described for these 
nuclei, even for cases with pronounced competing defor- 
mations, such as ^°Se and ™Se. The PCI dimensions are 
significantly lower the corresponding CI dimensions, as 
well as the corresponding computational effort. In addi- 
tion, in most cases, the low-lying projected basis states 
can provide some physical insight into the structure of 
the low-lying states. Finally, we calculated with the new 
method the low- lying 0+ states in ^^Ge and ^^Se that are 
relevant for the double beta decay of ''^Ge. The hope is 
that this method could be used some day to study the 
double beta decay of the strongly deformed ^^''Nd and 

238u 

Further improvements to the PCI method will include 
the extension of the formalism developed in Ref. [l^ to 
calculate electromagnetic transition probabilities. The 
new method uses different bases for different spins, which 
introduces additional complications. Other observables, 
such as spectroscopic amplitudes and DBD matrix ele- 
ments have to be worked out. Further improvement of 
the basis may be achieved for some cases that exhibit 
significant octupole deformation, which will require full 
projection on good parity. 
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